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ON  THE  LINEAR  STABILITY  OF  TWO-DIMENSIONAL  BARIUM  CLOUDS 

L  THE  INVISCID  CASE 

I.  Introduction 

The  study  of  barium  cloud  phenomenology  has  been  of  interest  to 
ionospheric  physicists  for  more  than  two  decades.  The  release  of  barium  in 
the  ionosphere,  followed  by  its  rapid  ionization,  leads  to  a  complex 
interaction  between  the  barium  plasma  and  the  background  ionospheric  plasma 
in  the  presence  of  ambient  electric  and  magnetic  fields.  One  of  the  most 
interesting  phenomena  observed  is  the  structuring  of  the  barium  cloud:  the 
"fingering"  or  "striating"  of  the  backside  of  the  cloud.  It  is  generally 
believed  that  this  structuring  is  produced  by  the  E  x  B  instability  [Linson 
and  Workman,  1970].  A  substantial  research  effort,  both  theoretical  and 
computational,  has  been  devoted  to  understanding  this  instability  and  its 
application  to  barium  cloud  structuring.  The  purpose  of  this  paper  is  to 
examine  the  linear  stability  of  the  steepened  backside  of  two-dimensional 
barium  clouds.  We  derive  a  relatively  simple  expression  for  the  growth 
rate  of  the  perturbations  on  these  backsides  which  are  functions  of  (1)  the 
parameter  M  which  Is  the  ratio  of  the  integrated  Pedersen  conductivity 
inside  the  barium  cloud  to  that  of  the  background  ionosphere;  (2)  the  shape 
of  the  cloud;  (3)  the  plasma  density  gradient  scale  length  L  of  the 
steepened  backside;  and  (4)  the  perturbation  wavenumber  k.  For 
sufficiently  small  kL,  the  growth  time  as  a  function  of  M  traces  out  the 
"U-shaped  curve"  observed  in  previous  studies  [Linson,  1975;  McDonald  et 
al. ,  1981;  Overman  and  Zabusky,  1980]. 

The  outline  for  the  remainder  of  this  paper  is  as  follows.  In  Section 
II  we  derive  the  growth  rate  of  a  two-dimensional  barium  cloud,  where  we 
approximate  the  steepened  backside  of  the  cloud  as  a  discontinuity.  In 
Section  III  we  relax  this  constraint  and  let  the  density  gradient  scale 
length  on  the  steepened  backside  be  finite.  In  Section  IV  we  relate  our 
results  to  the  "U-shaped  curve"  of  critical  Reynold's  number  seen  by 
McDonald  et  al.  [1981].  Finally,  in  Section  V  we  present  our  conclusions 
and  plans  for  future  work. 
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II.  Theory 


We  are  interested  in  deriving  the  approximate  linear  growth  rate  for 
perturbations  applied  to  the  steepened  backside  of  an  ionospheric  barium 
cloud.  To  achieve  this  we  need  several  earlier  results.  The  first  is  the 
equation  derived  by  Ossakow  and  Chaturvedi  [1978]  which  describes  the  rise 
velocities  of  plasma  depletions  in  the  shape  of  circular  or  ellipsoidal 
cylinders  in  the  presence  of  a  gravitational  field  perpendicular  to  a 
magnetic  field  3.  It  is  simple  matter  to  extend  their  equation  to  plasma 
enhancements  in  the  presence  of  an  electric  field  EL  (i.e.,  barium 
clouds).  The  second  result  we  shall  need  is  the  equation  for  the  linear 
growth  rate  of  the  gradient  drift  E  x  _B  instability  when  the  plasma  profile 
consists  of  a  single  discontinuity  in  plasma  density.  This  result  has  been 
derived  by  Huba  and  Zalesak  [1983].  Finally,  we  shall  need  to  know  the 
dependence  of  the  growth  rate  on  the  plasma  density  gradient  scale  length 
in  the  real-world  regime  between  the  two  extremes  of  a  density 
discontinuity  on  the  one  hand,  and  a  smooth  density  distribution,  where 
"local"  theory  should  hold,  on  the  other.  For  this  information  we  draw 
upon  the  work  of  Huba  et  al.  [1982]  and  a  suggestion  of  J.A.  Fedder 
[private  communication,  1983]. 

Consider  an  idealized  two-dimensional  barium  cloud  consisting  of  a 
uniform  density  inside  of  a  "waterbag"  of  ellipsoidal  shape,  acted  on  by  a 
uniform  external  electric  field  Eq  x  in  the  presence  of  a  magnetic  field  B 
z.  Let  the  ellipse  axis  in  the  direction  of  be  denoted  by  b  and  the 
other  axis  by  a.  Let  the  ion  density  of  the  barium  cloud  be  denoted  by 
and  that  of  the  background  ionosphere  be  denoted  by  nQ.  This  situation  is 
depicted  in  Fig.  1.  Then  the  results  of  Ossakow  and  Chaturvedi  [1978]  can 
be  extended  to  show  that  the  total  electric  field  E  inside  the  barium 
cloud  is  given  by 
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1  +  R 
1  +  MR 


(1) 


where 


r»  •  v  •  ^  \ 


M  i  (nb  +  no)/no  (3) 

Equation  (1)  is  the  first  result  we  shall  need  for  our  analysis. 

Now  consider  a  one-dimensional  barium  cloud,  infinite  and  invariant  in 
the  x  direction,  which  consists  of  a  single  jump  in  electron  density  from 
n^  +  nQ  to  nQ  at  y  »  0.  That  is,  ng  -  n^  +  nQ  for  y  <  0,  and  ng  *  nQ  for 
y  >  0,  where  ng  is  electron  density.  Again  we  impose  an  electric  field, 
denoted  here  by  ]£,  aligned  in  the  x  direction.  We  then  perturb  the 
interface  at  y  *  0  such  that  the  position  of  the  perturbed 

interface  n(x)  is  given  by  n  *  a  cos  kx  where  a  is  a  small  distance.  This 
situation  is  depicted  in  Fig.  2.  It  is  shown  in  Huba  and  Zalesak  [1983] 
that  this  interface  will  grow  in  amplitude  with  growth  rate  y  given  by 

Y  B  M  +  1 

Equation  (A)  was  derived  under  the  assumption  that  both  Hall  currents  and 
the  inertial  terms  in  the  ion  momentum  equation  are  negligible.  M  is  as 
defined  in  Eq.  (3).  An  alternative  derivation  of  Eq.  (A),  due  to  the  first 
author,  is  given  in  the  Appendix. 

Consider  a  steepened  barium  cloud,  as  depicted  in  Fig.  3.  If  we 
consider  only  the  tip  of  the  steepened  backside,  and  we  further  consider 
only  perturbations  of  sufficiently  long  wavelength  such  that  kL  is  small, 
where  k  is  the  perturbation  wavenumber  and  L  is  the  characteristic  scale 
length  of  the  density  gradient  on  the  backside  of  the  cloud,  then  the 
growth  rate  of  these  perturbations  should  be  well  approximated  by  Eq. 
(A).  However,  note  that  the  electric  field  £  seen  by  the  cloud  is  not  the 
externally  imposed  electric  field  E^,  but  rather  some  smaller  electric 
field  that  is  a  result  of  the  shielding  effects  of  the  two-dimensional 
cloud.  If  we  approximate  our  barium  cloud  as  an  ellipse  as  depicted  in 
Fig.  1,  then  Eq.  (1)  yields  precisely  thi3  smaller  electric  field: 

E  -  E  .  —  ~~  (5) 
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Substituting  this  field  into  Eq.  (4)  we  get  finally 

v  -  1.  CE°  <M  -  u  (R  •»■  n 

T  B  (M  +  1)  (MR  +  1) 

For  the  special  case  of  a  circular  shape  (R 

cE 


(6) 


1)  we  get 


2k 


_  M  -  1 

B  (M  +  l)2 


(7) 


Equation  (7)  traces  out  an  inverted  "U-shaped"  dependence  of  y  on 
which  we  relate  to  the  "U-shaped  curve"  of  McDonald  et  al.  [1981]  in 

Section  IV.  The  M  dependence  of  y  is  given  by  (M  -  1)/(M  +  l)2,  which 

maximizes  at  M  *  3.  A  plot  of  this  quantity,  as  well  as  its  two  factors, 

(M  -  1)/(M  +  1)  and  (M  +  1)  ^ ,  which  are  due  to  1-D  instability  and  2-D 
shielding  respectively,  is  shown  in  Fig.  4.  Note  that  the  falloff  in  y  for 
small  M  is  due  to  1-D  instability  effects  alone,  while  the  falloff  in  y  for 
large  M  is  due  to  2-D  shielding  alone.  Thus  we  conclude  that  the 


resistance  of  2-D  high  M  clouds  to  bifurcation  (see  McDonald  et  al. 


L  98 1 ) 


is  due  to  2-D  shielding. 

Equation  (7)  is  very  similar,  but  not  identical  to,  an  expression 
derived  by  Overman  and  Zabusky  [1980]  for  the  case  of  a  perfectly  circular 
"waterbag"  model  of  a  barium  cloud.  Their  linearized  analysis  showed  the 
results  of  perturbing  a  circular  cloud  of  radius  rQ  by  making  radial 
sinusoidal  displacements  of  its  boundary 


r  +  E  a  cos  m* 
o  m  T 

m 


(8) 


where  rQ  is  the  initial  radius  of  the  cloud,  r  is  the  perturbed  radius, 
ct^  is  the  perturbation  amplitude,  $  is  the  angle  measured  from  the  backside 
of  the  cloud,  and  m  is  an  integer  (mode  number).  They  found  that 

Note  that  if  we  write  our  perturbation  (8)  in  terms  of  arc  length  r  4> 
multiplied  by  a  wavenumber  k  then  cos  m6  ■  cos  k  rQ4>  and 


m/r  -  k  (10) 


4 


making  Eq.  (9)  yield  the  same  "growth  rate"  as  given  by  Eq.  (7),  except 

that  there  is  an  a  on  the  left  side  of  (9)  and  an  ct  ,,  on  the  right 

Ifl  TQ-rl 

side.  Thus  we  have  not  totally  recovered  the  results  of  Overman  and 
Zabusky  [1980]  for  the  circular  waterbag,  but  the  similarity  is  both 
striking  and  encouraging.  Furthermore,  our  formula  Eq.  (6)  exhibits 
several  advantages  over  the  Overman  and  Zabusky  [1980]  result.  For  one 
thing,  we  do  not  have  the  ambiguity  associated  with  using  the  term  "growth 
rate"  to  describe  the  rate  at  which  mode  m+1  supplies  energy  to  mode  m  (see 
Eq.  9).  Equation  (6)  is  a  genuine  growth  rate  in  the  usual  sense  of  the 
word.  Also,  precisely  because  we  have  derived  Eq.  (6)  as  the  product  of  a 

term  involving  global  effects,  i.e.,  the  overall  shape  of  the  cloud,  and  a 

term  involving  local  effects,  i.e.,  the  densities  and  electric  fields  at 
the  steepened  edge  of  the  cloud,  it  is  possible  to  generalize  the  results 
to  elliptical  clouds  (via  the  R  dependence)  and,  more  importantly,  to  the 
case  of  a  finite  density  gradient  on  the  backside  of  the  cloud,  as  we  shall 
see  in  Section  III. 

III.  Approximating  the  Growth  Rate  in  Cases  Where  the  Density  Gradient  is 

Finite 

Let  us  return  to  Eq.  (4).  This  equation  was  derived  under  the 
assumption  that  the  edge  of  the  barium  cloud  represented  a  true 
discontinuity  in  electron  density;  but  real  barium  clouds  consist  of  finite 
gradients.  How  accurately  does  Eq.  (4)  represent  the  case  of  a  sharp  but 
finite  density  gradient  and  how  are  we  to  treat  the  case  where  the  density 
gradient  is  not  even  sharp?  The  work  of  Huba  et  al.  [1983]  gives 


numerically  generated  growth  rates  for  a  wide  range  of  ratios  of 
perturbation  wavelength  to  gradient  scale  length.  However,  we  are 
interested  here  in  finding  a  simple,  closed  form  expression.  We  propose 
the  following,  after  a  suggestion  due  to  J.A.  Fedder  [private 
communication,  1983] 


,  cE  M  -  1  ,.T  .M+1 

K  b~  mtt  ’  L  *  inn: 


cE  T -1  ,  T  v  M  +  1 

y  -  —  L  ,  kL  > 
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where  L-^-  is  the  maximum  value  of  the  quantity  -j-j—  for  the  geometry 
shown  in  Fig.  2  wherein  we  have  replaced  the  discontinuity  by  a  finite 
transition  length.  Equation  (lib)  is  the  classical  result  for  the  gradient 
drift  instability  from  local  theory  (see  Llnson  and  Workman  [1970]).  The 
combined  Eq.  (11)  agrees  very  well  with  the  actual  numerical  results  shown 
in  Huba  et  al.  [1983]  when  allowance  is  made  for  the  fact  that  the  L  used 
in  that  paper  differs  from  the  one  used  here  by  a  factor  of  1.45.  A 
comparison  is  shown  in  Fig.  5. 

Combining  Eq.  (11)  with  Eq.  (1),  we  get  our  final  expressions  for  the 
growth  rate  of  instabilities  on  the  steepened  backside  of  a  barium  cloud, 
in  the  absence  of  diffusion: 


o  M  -  1  R  +  1 
B  M  +  1  MR  +  1 

:Eo  T-1  R  +  1 
B  MR  +  1 


M  +  1 
M  -  1 


M  +  1 
M  -  1 


(12a) 


(12b) 


An  obvious  caveat  to  the  above  is  that  the  perturbation  wavelength  X  = 
k/2n  must  be  smaller  than  the  total  diameter  of  the  barium  cloud  to  have 
the  concept  of  a  Fourier  mode  perturbation  to  make  any  sense  at  all. 

We  note  that  we  have  constructed  our  final  expressions  such  that  they 
not  only  very  closely  recover  the  results  of  Overman  and  Zabusky  [1980] 
(for  R  ■  1  and  L  *  0,  as  previously  discussed)  but  also  accurately 
reproduce  the  other  known  case  of  a  one-dimensional  plasma  distribution  (R 
■  0)  for  the  entire  range  of  kL. 


IV.  Relationship  to  the  U-Shaped  Curve 

Equation  (7),  when  viewed  as  an  expression  for  y  as  a  function  of  M, 
yields  a  curve  which  maximizes  at  M  ■  3  and  monotonically  decreases  away 
from  this  value  for  all  values  of  M  for  which  our  barium  cloud  represents  a 
true  enhancement  in  electron  density  (1  <  M  <  <*>) .  In  fact  the  growth  rate 
y  is  zero  for  M  ■  1  and  M  ■  «.  The  corresponding  curve  of  a  growth  time 
y  ^  versus  M  minimizes  at  M  *  3  and  goes  to  infinity  at  M  *  1  and  M  *  ®. 
Curves  of  y  *  versus  >  which  exhibit  this  qualitative  behavior  have  been 


dubbed  "U-shaped  curves"  by  several  researchers.  If  one  views  barium  cloud 
bifurcation  as  a  battle  between  the  gradient  drift  growth  rate  v  and 
diffusive  processes,  then  this  U-shaped  curve  can  be  related  to  a  similar 
U-shaped  curve  which  relates  the  amount  of  diffusion  D  needed  to  stop  the 
bifurcation  of  a  cloud  of  a  given  radius  rQ,  as  a  function  of  M.  Let  us 
assume  that  the  effect  of  diffusion  is  to  contribute  a  term  -  k  D  to  the 
growth  rate.  Let  us  further  approximate  our  barium  cloud  as  a  cylinder  (R 
=  1),  with  very  steep  edges  (L/rQ  <  l).  Then  our  growth  rate  is 
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M  -  1 

(M  +  l]2 


k2D 


(13) 


We  will  further  assume  that  bifurcation,  if  it  takes  place  at  all,  takes 
place  at  kL  ~  1.  The  bifurcation  condition  is  then  found  by  setting  y  =  0 
in  (13)  with  k  =  L-^: 
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Note  that  McDonald  et  al.  [1981]  got 
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so  the  functional  forms  are  the  same  but  our  Eq.  (16)  seems  to  state  that 
it  requires  150  times  more  diffusion  to  stop  bifurcation  then  was  observed 
by  McDonald  et  al.  [1981].  Whence  the  discrepancy?  Part  of  the  difference 
can  be  removed  by  noting  that  the  "finger"  of  McDonald  et  al.  [1981] 
corresponds  to  taking  R  *  in  which  case  we  would  get 


Pa 


M(M  +  1) 


For  large  M  this  removes  half  of  the  discrepancy.  Accounting  for  a  larger 
share  of  the  difference  is  ‘icDonald  et  al."s  [1981]  equation  for  computing 
L  (their  Eq.  (5)),  which  appears  to  give  L's  which  are  larger  than  the 
minimum  gradient  scale  length  we  have  assumed  here.  However,  we  do  not 
believe  this  effect  could  account  for  a  factor  greater  than  ten,  which 

still  leaves  a  discrepancy  of  more  than  seven.  That  is,  the  model  we 

propose  here  predicts  a  factor  of  at  least  seven  greater  diffusion 

necessary  to  stop  bifurcation  than  is  given  by  McDonald  et  al.  [1981].  The 
U-shaped  dependence  on  M,  however,  is  the  same.  We  hope  to  be  able  to 
resolve  these  quantitative  differences  in  the  sequel  to  this  paper 

[Zalesak,  1983],  where  the  full  effects  of  diffusion  are  considered. 

V.  Conclusions 

We  have  derived  an  expression,  Eq .  12,  that  we  believe  to  be  a  good 
approximation  to  the  actual  growth  rate  of  perturbations  applied  to  the 
steepened  backside  of  a  two-dimensional  ionospheric  barium  cloud,  for  the 
case  of  no  electron  or  ion  diffusion.  Diffusion  plays  a  large  role  in  the 
behavior  of  two-dimensional  barium  clouds,  however,  and  affects  the  growth 
rate  in  two  ways: 

(i)  diffusion  will  determine  the  maximum  density  gradient  that  can 
exist  at  the  backside  of  the  barium  cloud; 

(ii)  diffusion  will  provide  a  damping  mechanism  for  any  perturbation 
which  attempts  to  grow  on  the  backside  of  the  cloud.  We  have 
already  treated  this  effect,  at  least  approximately,  in  Section 
IV. 

These  two  combined  effects  of  diffusion  are  treated  in  the  sequel  to  this 
paper  [Zalesak,  1983]. 


Appendix.  Derivation  of  the  E  x  B  Instability  Growth  Rate  for  a 
Discontinuous  Interface 


The  geometry  we  shall  consider  is  depicted  in  Fig.  6.  The 


differential 

equations  describing  this  system  are  [McDonald  et  al., 

3n 

1981]  : 

JT  *  '  0 

(Al) 

V*[neV(j)J  =  E*Vng 

(A2) 

where 

v  =  -  ~  V<j>  x  B 

B2 

(A3) 

Here  ng  is  the  electron  density,  which  is  n'  for  y  >  n(x)  and  n  for 
y  <  p(x),  c  is  the  speed  of  light,  j}  is  the  magnetic  field,  E  is  the 
imposed  electric  field,  and  <j>  is  the  induced  electrostatic  potential.  Note 
from  Eq.  (A3)  that  we  have  placed  ourselves  in  a  reference  frame  drifting 
with  the  £  x  _B  velocity.  Our  two-dimensional  plasma  resides  in  the  x-y 

A 

Cartesian  plane.  B  is  aligned  along  z. 


If  E  *  E  x,  and  we  ignore  Hall  currents,  then  the  geometry  depicted  in 
Fig.  6  is  an  equilibrium  configuration  with  n  ”  0.  We  now  apply  a 
perturbation  to  the  position  of  the  interface,  denoted  by  n(x),  of  the  form 


! 

* 


a  cos  kx  e 


iat 


(A4) 


.  and  ask  what  the  response  of  the  electrostatic  potentials  in  the  two 
I  regions  will  be.  Noting  that  Vn  *  Vn'  =  0,  since  the  fluid  is 
j  incompressible  and  we  have  only  moved  the  interface,  Eq.  (A2)  becomes 

i 


V2f  -  0  ,  y  >  n(x) 

V2*  -  o  ,  y  <  n(x) 


If  we  assume  solutions  of  the  form 


<b  m  G'(y)  sin  kx  e 


iat 


y  >  n(x) 


(A5) 

(A6) 


(A7) 
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6  *  G(y)  sin  kx  e 


y  <  r(x) 


(A8) 


(Note  we  are  assuming  a  specified  phase  relationship  between  n  and  *,  which 
will  later  be  shown  to  satisfy  our  boundary  conditions),  then  demanding 


f  vanish  at  y  *  +  ®  and  $  vanish 
<j>'  *  C'  e  ^  sin  kx  e^0t , 
d>  *  C  e"*^  sin  kx  e*0t, 


y  »  -  *  yields 


y  >  n(x) 

(A9) 

y  <  n(x) 

(A10) 

We  shall  later  see  that  a  is  pure  imaginary,  meaning  that  the  perturbation 
does  not  propagate  in  the  x  direction.  Therefore  the  movement  of  the 
interface  is  determined  by  the  y-components  of  the  velocity  fields  at  the 
interface,  which  must  be  the  same  for  the  interface  to  be  well  defined  (no 
interpenetration  or  cavitation): 


but 


v  at  y 

y 


n 


v 

y 


c  3_4 
B  3x 


v  * 

y 


c  3i>' 
B  9x 


(All) 


so 

B  3n  ,  _ .  -ky  ,  iat  i 
—  -rr  =  kC  e  cos  kx  e 
c  3t  y  -  n 

*  kC  e"*^^  cos  kx  e*0t| 

y  -  n 

Since  n  is  assumed  to  be  arbitrarily  small  we  have 

c-  -  C 


(A12) 


(A13) 


but  by  Eq.  (A4) 


3n 


ia  a  cos  kx  e 


iat 


(A14) 


so  comparing  (A12)  and  (A14)  we  gee 

i  |  a  a  -  kC'  -  kC  (A15) 

Now  we  have  one  last  constraint  to  satisfy,  namely  that  7*J  -  0.  We 
will  use  the  integral  form  of  this  constraint,  that  is 


Jl 


S(V)-‘n  dS 


(A16) 


for  any  control  volume  V  bounded  by  surface  S(V).  Equation  (A16)  obviously 
holds  for  any  control  volume  not  containing  the  interface  n(y),  since  n  *  a 
constant  and  7^$  *  0.  Consider  our  interface  and  the  control  volume 
depicted  in  Fig.  7.  J^,  Jg,  Jg  and  J^,  are  the  total  currents  flowing 
through  the  left,  right,  bottom  and  top  surfaces  respectively  of  the 
control  volume.  We  define  an<*  ’"'l  s  n(  x^J  . 


where 


JJ  i*n  ds  -  J  -  J  +  J  -  J  -  0 


J  -  -  n 
T  3y 


“kyT  XT? 

+  n'kCe  I  sin  kx  dx 

*  V 


a*  kyR  xr 

j  .  -  n  -  -  nkCe  f  sin  kx  dx 
B  3y  JxL 

JR  •  JR  +  JR 


TII  #  R  f_ 
JR  5  iy  ntE 
yB 

jjf  5  J71  n'lE 
R  \  1 


36,,. 

3l^dy 


(A17) 

(A18) 

(A19) 

(A20) 

(A21) 

(A22) 


3  <b 

Since  both  - —  and  the  integration  path  length  [y  -  n„)  are  perturbation 

dX  1  K' 

quantities  in  Eq.  (A22),  the  integral  involving  n"  36' /3x  will  be  quadratic 
in  the  perturbation  and  hence  negligible.  We  then  get 


J 


I 

R 


yT 

I  n'E  dy 

% 


(A23) 
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Similarly 


IT 

J„  *  |  a  E  dy  *  h  -  yJnE 

R  Jy.  3  y  R  'B- 

D 

JR  ■  JR  +  J”*  V'E  '  VE  +  V”  -  >E 

Using  similar  reasoning  we  get 

JL  »  yTn'E  -  yBn£  +  njn  -  n')E 
JR  -  JL  =  (n  -  n  )E(nR  -  nL) 

We  now  let 


(A24) 

(A25) 


(A26) 

(A27) 


(*R  “  \)  *  0  (A28) 

Then 

JR  -  JL  -  -  n')E  |1  (>R  -  xj  (A29) 

but  from  Eq.  (A4) 

|j  -  -  ctk  sin  kx  eiat  (A30) 

so 

JR  -  *  -  (n  -  n')E  a  k  sin  kx  eiat[xR  -  x^j .  (A31) 

Using  Eq.  (A28)  and  noting  that  y^,  and  y^  are  perturbation  quantities  we 
can  set 

XR 

/  sin  kx  dx  *  sin  kx(xD  -  xT J  (A32) 

XL  K  L 

-ky  +ky 

e  -  e  =  1  (.433) 


(  A34) 


and  hence  Eqs.  (A18)  and  (A19)  become 

JT  *  +  n'kC  sin  kx(xR  -  x^J 


J_  =  -  n  kC  sin  kx(x„  -  xT ) 
B  R 


JT  ”  JB  *  +  kC  Sin  kx^xn  "  xt) 


R  V 


(A35) 

(A36) 


Combining  Eq.  (A36),  (A31),  and  (A17)  and  canceling  the  common  term 
k  sin  kx  elcrt(xR  -  x^J  we  get 


(n  +  n)C  -  (n  -  n')E  a  *  0 


(A37) 


Now  substituting  Eq.  (A15)  into  (A37) 


(n  +  n)i 


Baa 

KT 


(n  -  n  )E« 


.  n  -  n  ,  cE 

a  =*  -  i  — - -  k  — 

n  +  n  B 


n  -  n  ,  cE 

y  ■ - - - r  k  —— 

n  +  n  B 


Defining  M  =  n/n" 

Y  -  1  k  CE 

Y  M  +  1  B 

which  corresponds  to  the  result  obtained  by  Huba  and  Zalesak  [1983]. 


(A38) 

(A39) 


(A40) 


£  . 
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.  Idealized  "waterbag”  elliptical  barium  cloud  geometry.  The  total 
electron  density,  barium  cloud  ion  density,  and  background 
ionosphere  ion  density  are  denoted  by  ng,  n^,  and  nQ 
respectively.  _B  is  the  ambient  magnetic  field,  and  E^  is  the 
externally  imposed  electric  field.  E  is  the  constant  shielded 
electric  field  inside  the  cloud.  a  and  b  are  the  major  and  minor 
axes  of  the  ellipse  respectively. 


X 


Geometry  of  1-D  barium  cloud  configuration  for  stability  analysis 
of  discontinuous  interface. 


Isodensity  contours  of  a  steepened  2-D  barium  cloud.  The  box 
indicates  that  region  of  space  over  which  we  will  apply  our  1-D 
stability  analysis  by  approximating  the  density  gradient  there  as 
a  1-D,  slowly  evolving  structure. 


- D  GROWTH  RATE  ~ 


TOTAL  GROWTH  RATE  ~ 


(M  +  l)‘ 


2  -  D  SHIELDING  ~  (M  +l)‘ 


Plot  of  the  M-dependence  of  the  growth  rate  given  by  Eq.  (7),  M  - 
1/ (M  +  1)  .  Also  shown  are  its  two  component  parts,  (M  -  1)/(M  + 
1)  and  (M  +  1)“*,  due  to  1-D  instability  and  2-D  shielding 
results  respectively. 


w 


Plot  of  normalized  growth  rate  versus  kL  for  the  case  of  the  M  = 
39  cloud  considered  by  Huba  et  al.  [1982].  The  solid  line 
depicts  the  simple  model  shown  here,  while  the  circles  represent 
numerical  results  from  Huba  et  al.  [1982]. 


A 

E  =E  x 


FLUID  INTERFACE  y  *  ^  lx)  =  a  cos  kx 


n',*' 


n  ,  4> 


-y,-0 


•CONTROL  VOLUME  V 


•  B  s  B  z 


Geometry  of  1-D  barium  cloud  configuration  for  stability 
analysis,  as  in  Fig.  2.  n  and  6  are  the  plasma  density  and 
electrostatic  potential,  respectively.  Primed  and  unprimed 
variables  denote  quantities  above  and  below  the  interface 
y  ■  n(x),  respectively.  Shown  is  the  control  volume  used  to 
evaluate  V»J. 
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01CY  ATTN  OPR  ALVA  T.  STAIR 

01CY  ATTN  PHD  JURGEN  BUCHAU 

01CY  ATTN  PHD  JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
01CY  ATTN  SUL 

01CY  ATTN  CA  ARTHUR  H.  GUENTHER 
01CY  ATTN  NTYCE  1LT.  G.  KRAJEI 

AFTAC 

PATRICK  AFB,  FL  32925 
01CY  ATTN  TF/MAJ  WILEY 
01CY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  AAD  WADE  HUNT 
01CY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
01CY  ATTN  AFRDQ 


01CY  ATTN  CODE  DF-14  R.  BUTLER  HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
OFFICER  OF  NAVAL  RESEARCH  DEPARTMENT  OF  THE  AIR  FORCE 

ARLINGTON,  VA  22217  HANSCOM  AFB,  MA  01731 

01CY  ATTN  CODE  465  01CY  ATTN  J.  DEAS 

01CY  ATTN  CODE  461 

01CY  ATTN  CODE  402  HEADQUARTERS 

01CY  ATTN  CODE  420  ELECTRONIC  SYSTEMS  DI VISION/YSEA 

01CY  ATTN  CODE  421  DEPARTMENT  OF  THE  AIR  FORCE 

HANSCOM  AFB,  MA  01732 
01CY  ATTN  YSEA 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COKMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  ADWATE  MAJ  3 RUCK  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL 

COMMANDER 


DEPARTMENT  OF  ENERGY 
AL3UQ’JERQUE  OPERATIONS  OFFICE 
P.O.  30X  5400 
ALBUQUERQUE,  NM  37115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG4G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  308 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


OICY 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

ZINN 

OICY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTER VEL 

OICY 

ATTN 

D.  : 

SAPPENFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
AL3UQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 


ROME  AIR  DEVELOPMENT  CENTER,  AFSC  SANDIA  LABORATORIES 

HANSCOM  AFB,  MA  01731  LIVERMORE  LABORATORY 

OICY  ATTN  EEP  A.  LORENTZEN  P.O.  BOX  969 

LIVERMORE,  CA  94550 

DEPARTMENT  OF  ENERGY  OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 

LIBRARY  ROOM  C-042  OICY  ATTN  DOC  CON  FOR  T.  COOK 

WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ  OFFICE  OF  MILITARY  APPLICATION 

DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  COS  DR.  YO  SONG 


OTHER  GOVERNMENT 

DEPARTMENT  OP  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

01CY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  4  INFO  ADMIN 
BOULDER,  CO  80303 
01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
OICY  ATTN  L.  SLOAN 
OICY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HKES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  HILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  HA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  HA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  6  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CT  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DAS I AC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINXA3IT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-1 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  COR? 

ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  PL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
01CT  ATTN  J.  MOULE 
01CT  ATTN  GEORGE  MROZ 
01CT  ATTN  W.  OLSON 
01CT  ATTN  R.W.  HALPRIN 
01CT  ATTN  TECHNICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  R.  BOGUSCH 
OICY  ATTN  R.  HENDRICK 
OICY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FA  JEN 
OICY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGMIRE 
OICY  ATTN  B.  WHITE 


BEDFORD,  HA 
OICY  ATT 
OICY  ATT 


HA  01730 

ATTN  JOHN  MORGANSTERN 
ATTN  G.  HARDING 


ATTN  C.E.  CALLAHAN 


MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  Vk  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 


PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 


PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


MISSION  RESEARCH  CORP. 

OICY 

ATTN 

FORREST  GILMORE 

1720  RANDOLPH  ROAD,  S.E. 

OICY 

ATTN 

WILLIAM  B.  WRIGHT, 

ALBUQUERQUE,  NEW  MEXICO  87106 

OICY 

ATTN 

ROBERT  F.  LELEVIER 

OICY  R.  STELLINGWERF 

OICY 

ATTN 

WILLIAM  J.  KARZAS 

OICY  M.  ALME 

OICY 

ATTN 

H.  ORY 

OICY  L.  WRIGHT 

OICY 

ATTN 

C.  MACDONALD 

OICY 

ATTN 

R.  TURCO 

MITRE  CORPORATION,  THE 

OICY 

ATTN 

L.  DeRAND 

P.O.  BOX  208 

OICY 

ATTN 

W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDR0Z1AN 


RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 


RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 


IvV V  jvV.V. ' 
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SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


OICY 

ATTN 

DONALD  NEILSCN 

OICY 

ATTN 

ALAN  BURNS 

OICY 

ATTN 

G.  SMITH 

OICY 

ATTN 

R.  TSUNODA 

OICY 

ATTN 

DAVID  A.  JOHNSON 

OICY 

ATTN 

WALTER  G.  CHESNUT 

OICY 

ATTN 

CHARLES  L.  RINO 

OICY 

ATTN 

WALTER  JAYE 

OICY 

ATTN 

J.  VICKREY 

OICY 

ATTN 

RAY  L.  LEADABRAND 

OICY 

ATTN 

G.  CARPENTER 

OICY 

ATTN 

G.  PRICE 

OICY 

ATTN 

R.  LIVINGSTON 

OICY 

ATTN 

V.  GONZALES 

OICY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 
OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  4  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLESUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D /  STOCKKELL 
SNTF/1575 

VT3IDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

Dr.  P.  MANGE  -  CODE  4101 
Dr.  P.  COODMAN  -  CODE  4130 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14350 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUB IN 

DR.  N.  MAYNARD  -  CODE  696 
COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZ'JBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 

ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  CLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NWI6-262 
CAMBRIDGE,  MA  02139 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  P.0KR3AUGH 
DR.  L. A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  E.  FRIEMAN 


yj 


STAFFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY  . 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AX  99701 
DR.  L.E.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  BARASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 
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